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ABSTRACT 

Recent studies suggest that when magnetohydrodynamic (MHD) turbulence is excited by stirring a plasma 
at large scales, the cascade of energy from large to small scales is anisotropic, in the sense that small-scale 
fluctuations satisfy the inequality k\\ <C k±, where £|| and k± are, respectively, the components of a fluctuation's 
wave vector || and _L to the background magnetic field. Such anisotropic fluctuations are very inefficient at 
scattering cosmic rays. Results based on the quasilinear approximation for scattering of cosmic rays by anisotropic 
MHD turbulence are presented and explained. The important role played by molecular-cloud magnetic mirrors in 
confining and isotropizing cosmic rays when scattering is weak is also discussed. 



1. INTRODUCTION 

In diffusion models of Galactic cosmic -ray propagation, cos- 
mic rays are scattered by small-scale fluctuations in the in- 
terstellar magnetic field. For cosmic -ray energies below ~ 
10 2 — 10 3 GeV, these small-scale fluctuations can arise from 
resonant waves that the cosmic rays generate themselves. At 
higher energies, it is believed that self-confinement is not possi- 
ble, because the growth rates of the resonant modes become too 
small in comparison to the rates at which the modes are damped 
(Cesarsky 1980, Berezinskii et al. 1990). For cosmic rays with 
energies above 10 2 — 10 3 GeV, scattering can result from turbu- 
lence that is generated by large-scale stirring of the interstellar 
medium (ISM), which results in a cascade of magnetic energy 
from large to small scales. Recent studies of magnetohydro- 
dynamic (MHD) turbulence, however, find that the small-scale 
fluctuations resulting from a turbulent cascade satisfy the in- 
equality <C k±, where k\\ and k± are, respectively, the compo- 
nents of a fluctuation's wave vector || and _L to the background 
magnetic field Bo (e.g., Goldreich & Sridhar 1995). This paper 
explains why fluctuations with this anisotropy are very ineffi- 
cient at scattering cosmic rays, and presents results from quasi- 
linear theory for the scattering mean free paths resulting from 
anisotropic MHD turbulence. The role of molecular clouds in 
confining and isotropizing cosmic rays when scattering is weak 
is also discussed. 

2. ANISOTROPIC MHD TURBULENCE AND THE 
GOLDREICH- SRIDHAR SPECTRUM 

Early studies of MHD turbulence assumed that when a 
plasma is stirred on some large scale I, the cascade of energy 
from large scales to small scales proceeds isotropically in k- 
space (Kraichnan 1965). More recent studies, however, find 
that energy cascades efficiently to large values of k±, but not 
very efficiently to large values of k\\ (e.g., Goldreich & Sridhar 
1995, Shebalin et al. 1983). (Note: if the mean-magnetic field 
is weaker than the fluctuating magnetic field, then within any 
stirring-scale cell of volume Z 3 there is a preferential field direc- 
tion which can be thought of as the background field Bo for all 
of the small-scale fluctuations within that cell. Local anisotropy 
is then determined relative to the direction of Bo within each 
stirring-scale cell.) 



Goldreich & Sridhar (1995) have proposed an inertial-range 
power spectrum for strong anisotropic MHD turbulence in 
which the magnetic field fluctuations are comparable to Bo: 

E B (k ± ,^ki^r^ g (^^, (i) 

where the dimensionless function g(x) is ~ 1 for |x| < 1 and 
rapidly approaches for |x| ^> 1. Evidence in support of equa- 
tion (1) has been found in direct numerical simulations of MHD 
turbulence (Maron 2000, Cho & Vishniac 2000). In this spec- 
trum, there is only power at small scales when k\\ < £"^ 3 /~ 1//3 . 
In this region of k-space, in which turbulence is excited, the 
linear incompressible Alfven-wave period (&||Va) is greater 
than the nonlinear energy-transfer time (k^v^ 1 , where va = 
Bo/^/4np is the Alfven speed, p is the mass density of the 
medium, and ~ VA{k±l)~ 1 ^ is the rms velocity fluctuation 
on a perpendicular scale of kj 1 . Because of this inequality, the 
Iroshnikov-Kraichnan mechanism for slowing energy-transfer 
does not apply (Kraichnan 1965). 

3. WHY SCATTERING IS WEAK IN ANISOTROPIC 
TURBULENCE 

If the magnetic power spectrum is isotropic in k-space or 
possesses slab symmetry (in which the wave vectors of fluc- 
tuations are || to Bo), then cosmic-ray scattering is dominated 
by magnetostatic gyroresonant interactions, in which the cos- 
mic ray and fluctuation satisfy the resonance relation 

fc||V||=nfl, (2) 

where V|| is the component of a cosmic ray's velocity along Bo, 
n is a non-zero integer (n = ±1 for slab symmetry), and Q. is 
the cosmic ray's gyrofrequency. In the general resonance rela- 
tion, the linear wave frequency co appears on the left-hand side 
of equation (2). It is neglected here since for incompressible 
Alfven waves co = k\\\>A, which is <C fc||V|| unless the angle be- 
tween a cosmic-ray's velocity vector v and Bo, the pitch angle, 
is very close to 90°. 

If a cosmic ray's pitch angle isn't too close to 0°, 90°, or 
180°, then equation (2) implies that the fluctuations that dom- 
inate scattering satisfy k\\ ~ p -1 , where p = Vj_/i2 is a cosmic 
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ray's gyroradius and v± is the component of v _L to Bo. How- 
ever, if p 1 » where I is the scale at which the turbulence 
is stirred, and if the turbulence has the type of anisotropy de- 
scribed by equation (1), then the only fluctuations with k\\ <~ p _1 

have k ± » p _1 , as depicted in figure 1. But if k ± p > 1, then 
during a single gyro orbit a cosmic ray traverses many uncor- 
rected turbulent fluctuations of the required k^. The contribu- 
tions from these different fluctuations tend to cancel, resulting 
in highly inefficient scattering. 

Scattering rates have been calculated for turbulence with a 
power spectrum described by equation (1) in the quasilinear ap- 
proximation (Chandran 2000a). If 
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then when e 3 / 2 <C (— lne)8 the coefficient of spatial diffusion 
along the magnetic field resulting from the quasilinear scatter- 
ing rates is given by (Chandran 2000a) 

(5) 



This value is far too large to explain confinement of cosmic rays 
to the Galaxy. Thus, some mechanism besides turbulence de- 
scribed by equation (1) must be invoked to explain cosmic-ray 
confinement. For cosmic rays with energies less than 10 2 -10 3 
GeV, waves excited by the cosmic rays may provide the con- 
finement. For higher energy cosmic rays for which self con- 
finement does not appear possible (Cesarsky 1980, Berezinskii 
et al. 1990), molecular-cloud magnetic mirrors may play an 
important role. It should be noted that equation (5) applies to 
both cosmic -ray nuclei and electrons. 

4. WHY MOLECULAR CLOUDS CAN HELP CONFINE COSMIC 
RAYS WHEN SCATTERING IS WEAK 

Molecular-clouds are characterized by a range of sizes, 
masses, and densities. Cloud mass spectra obey power-law 
scalings over several decades of cloud masses, and a power-law 
mass-size relation also holds over a range of scales (Elmegreen 
& Falgarone 1996, Heithausen et al. 1998, Blitz & Williams 
1997). Elmegreen (1997) has proposed a useful quasi-fractal 
model for this hierarchy of structure, with smaller, denser ob- 
jects nested within structures that are larger and more dif- 
fuse. Using straightforward rules to generate fractal struc- 
tures, Elmegreen estimates through numerical modeling that a 
line-of-sight through a fractal molecular cloud complex has a 
50 ± 10% chance of entering dense molecular material, and a 
50 ± 10% chance of passing through a hole in the fractal com- 
plex filled with diffuse matter. Elmegreen interprets the stan- 
dard 8 large absorption lines per kpc (Blaaw 1952) as evidence 
for an average of 3 fractal cloud complexes per kpc along a typ- 
ical line of sight. Because of the 50% see-through probability, 
however, photons can travel <~ 600 pc without entering molecu- 
lar material. The 8 absorption lines per kpc are then not spaced 
at even distances along a line of sight, but rather are clustered 
in groups of <~ 5 per cloud complex. 

Because magnetic field lines are focused into strong-field re- 
gions, and because the magnetic field is stronger within molec- 
ular clouds than in the ICM, a magnetic field line passing 
through a cloud complex has a higher probability of entering 
molecular material than a straight line of sight (figure 2). If the 



mean field strength within the molecular material is m times 
the typical field strength in the ICM, and if 50% of the lines- 
of-sight through the complex intersect molecular material, then 
one would expect a fraction P on the order of m/(m + 1) of the 
magnetic flux (and, therefore, field lines) through the complex 
to pass through molecular material. The line-of-sight average of 
a cloud's magnetic field can be obtained through Zeeman split- 
ting. Troland & Heiles (1996) report Zeeman measurements of 
molecular-cloud field strengths ranging from 9 /jG to 120 /jG. 
The typical field strength in the ICM is ^4 — 5 /jG (Zweibel 
& Heiles 1997). The ratio m is thus fairly large, and P is close 
to 1 . Moreover, due to magnetic focusing, the spacing of cloud 
complexes along field lines may be smaller than the <~ 300-pc 
spacing of complexes along straight lines. On the other hand, 
tangling of field lines in the ICM may tend to increase the spac- 
ing of successive cloud complexes along a given field line as 
measured along that field line. It seems reasonable, however, to 
take the typical distance between cloud complexes as measured 
along a field line to be 

Wrcioud — 300 pc. (6) 

Once inside a molecular cloud, field lines tend to be focused 
into dense clumps where the field strength is larger than the av- 
erage field strength in the cloud. A range of observations sug- 
gests that at particle densities n above 10 2 cm~ 3 , the magnetic 
field strength/? scales as (Vallee 1997) 



(7) 



According to Heithausen et al. (1998) 
Mccr 2 ' 3 , 



n cx r 7 , and 



(8) 
(9) 

dN/dr^r 30 , (10) 

where M is the mass of a clump of linear dimension r, and 
(dN /dr)Ar is the number of clumps with linear dimension r in 
the interval (r,r+ Ar). Since the flux through a clump is ~ Br 2 , 
equations (7) through (10) imply that the total flux 4> through 
clumps with linear dimension between r and 2r scales as 

<t>oc r -° 35 . (11) 

That is, there is more flux through the smaller, denser clumps 
than through the cloud complex as a whole. This means that 
as a single field line passes through a complex, it must on aver- 
age pass through several of the densest clumps described by 
the power-law scalings, clumps in which the magnetic field 
strength is large. 

If Bicm is the field strength in the ICM and B max is the maxi- 
mum field strength encountered by a cosmic ray in a molecular- 
cloud complex, then the cosmic ray will be magnetically re- 
flected by the cloud complex provided that its pitch-angle co- 
sine ^icm = V||/v in the ICM as it approaches the complex sat- 
isfies the inequality 



ISicmI < yJl-X' 1 , 



where 



B„ 



Bicm 



(12) 



(13) 



The above discussion suggests that % ^> 1, and thus molecular 
clouds can magnetically reflect a large fraction of cosmic rays. 

When the scattering mean-free path Ky /v is much greater 
than /intercioud. cosmic rays travel between molecular clouds 
without significant scattering. Phase-space orbits of cosmic 
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Fig. 3.— 



rays in this weak-scattering limit are depicted in figure 3, where 
£ = V||/v. Particles with |^| close to 1 are not magnetically re- 
flected and can pass through molecular clouds. Trapped cosmic 
rays move on closed orbits in the x-t, plane, where x is distance 
along a field line. 

The way in which molecular-cloud magnetic mirrors af- 
fect cosmic -ray transport depends upon cosmic-ray energy and 
also the efficiency of cosmic -ray scattering. Approximate 
coefficients of diffusion perpendicular to the Galactic plane 
have been calculated for several different propagation regimes 
(Chandran 2000b). Broadly speaking, cosmic rays can escape 
from magnetic traps in one of two ways, by scattering into the 
passing region of phase space and then traveling along the mag- 
netic field through a molecular cloud, or by drifting perpendic- 
ular to the magnetic field. The second case is depicted in figure 
4: a cosmic ray initially trapped between clouds A and B can 
drift perpendicular to the magnetic field and end up trapped be- 
tween clouds A and D. 

Because molecular clouds are confined to the Galactic disk, 
it is possible that they have no affect on cosmic-ray propaga- 
tion in the halo. On the other hand, if the magnetic field lines 
in the halo possess numerous arcs that are anchored down to 
the disk on either end (analogous to closed field lines in the so- 
lar corona), then molecular cloud magnetic mirrors may affect 
propagation in the halo to some extent. 

5. ISOTROPIZATION OF COSMIC RAYS BY 
MOLECULAR-CLOUD MAGNETIC MIRRORS 

One of the sources of cosmic-ray anisotropy is the flow of 
cosmic rays along the magnetic field. Because molecular-cloud 
magnetic mirrors impede this flow, they reduce the level of 
anisotropy for any given level of weak scattering. Values of the 
harmonics of the cosmic -ray distribution function (as functions 
of the scattering rate, %, and /i nterc ioud) are given by Chandran 
(2000b) under the assumption that the pitch-angle scattering 
frequency is weak and independent of pitch angle. 



6. IMPLICATIONS OF MOLECULAR-CLOUD MAGNETIC 
MIRRORS FOR DIFFUSE GAMMA RADIATION AND 
SECONDARY PRODUCTS 

When scattering is weak, the density of cosmic rays within 
molecular clouds n c ioud is determined by two competing ef- 
fects. On the one hand, cosmic rays are reflected as they ap- 
proach cloud complexes, which tends to reduce n c ioud- On the 
other hand, magnetic field lines are brought closer together in 
high-field regions, which acts to increase n c i ou d since cosmic 
rays travel primarily along the magnetic field. It can be shown 
that when energy losses are neglected, these two effects can- 
cel (Chandran 2000b). This point is important since if n c ioud 
were in fact less than «icm, there would be a corresponding re- 
duction in spallation and diffuse gamma radiation for a fixed 
average energy density of cosmic rays throughout the Galaxy. 
(For sufficient ionization losses at low cosmic-ray energies, it 
should be noted that the value of « c ioud can be reduced below 
the cosmic -ray density in the intercloud medium «icm-) 

The two competing effects described above are illustrated 
graphically in figure 5. Each of the two narrow flux tubes in 
figure 5 has a bounding surface that is everywhere parallel to 
the magnetic field. The cross-sectional area of each tube is pro- 
portional to (l/B), where B is the field strength. Since mo- 
tion perpendicular to the magnetic field is suppressed, the cos- 
mic rays within each flux tube to a good approximation remain 
within their respective flux tubes as they move along the field. 
Because of magnetic mirroring, the number of cosmic rays per 
unit length within a flux tube decreases in high-field regions. 
However, because the cross-sectional area of the flux tube also 
decreases, the number of cosmic rays per unit volume stays the 
same. 

7. DOES THE MODEL FIT THE DATA? 

At this stage, it is difficult to determine from observations 
whether molecular clouds play a role in cosmic ray confine- 
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Fig. 5.— 



ment. Although the observed energy dependence of the cosmic- 
ray path length at cosmic ray energies < 10 2 GeV provides 
important information on propagation at energies < 10 2 GeV, 
almost nothing is known about the path length at the ener- 
gies above 10 — 10 3 GeV at which self-confinement appears 
to break down and at which confinement may depend upon 
molecular clouds. There appear to be three main possibilities. 
First, molecular clouds may help confine cosmic rays at ener- 
gies above 10 2 — 10 3 GeV as described in this paper. Second, 
the arguments that self-confinement breaks down at energies 
above 10 2 — 10 3 GeV may be incorrect. Third, interstellar tur- 
bulence generated by large-scale stirring may possess features 
not described by the Goldreich-Sridhar theory that allow for 
stronger scattering. 

8. CONCLUSION 



Recent investigations into MHD turbulence are providing 
new and important results on the anisotropy of the small-scale 
fluctuations that result from a cascade of magnetic energy from 
large to small scales. As discussed in this paper, anisotropic 
small-scale fluctuations are inefficient at scattering cosmic rays. 
For cosmic rays with energies less than 10 2 10 3 GeV, resonant 
waves excited by streaming cosmic rays are believed to be suf- 
ficient to confine cosmic rays to the Galaxy regardless of the 
nature of the cascade in MHD turbulence. At higher energies, 
however, it is believed that such self-generated waves are insuf- 
ficient. Thus, if scattering by the turbulence that is generated 
by large-scale stirring of the ISM is inefficient, then some ad- 
ditional mechanism is needed to confine and isotropize cosmic 
rays at energies above ~ 10 2 — 10 3 GeV. Such a mechanism 
may be provided by molecular-cloud magnetic mirrors. 
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